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The nonlinear dynamics of a Bose-Einstein condensate (BEC) of dipolar excitons trapped in an 
external confining potential in coupled quantum wells is analyzed. It is demonstrated that under 
typical experimental conditions the dipolar excitons BEC can be described by a generalized Gross- 
Pit aevskii equation with the local interaction between the excitons, which depends on the exciton 
distribution function. It is shown that, if the system is pumped at sufficiently high frequencies, a 
steady turbulent state can be formed. 
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In the last decade, the nonlinear dynamics of the ex- 
citations in semiconductor heterostructures coupled with 
laser radiation attract attention because of promising po- 
tential applications in electronics and photonics, includ- 
ing, e.g., design of thresholdless lasers, optical comput- 
ing and quantum computing One of the exam- 
ples of the systems where the excitations demonstrate 
essentially quantum behavior is a Bose-Einstein conden- 
sate (BEC) of excitons in semiconductors ^^-tSi]- The dy- 
namics of dipolar excitons formed by spatially separated 
charges in coupled quantum wells (CQWs) in semicon- 
ductor heterostructures at helium temperatures attract 
attention due to the relatively long exciton lifetime com- 
pared to excitons in a single quantum well Re- 
cently, it was predicted that dipolar excitons can form a 
superfluid in a two-layer graphene 0. 

In this paper we focus on non-stationary dynamics in 
a dipolar exciton BEC in QCWs in an in-plane trapping 
potential. The trapping potential is essential for the con- 
densate formation at finite temperatures and it can be 
formed in GaAs structures by applying inhomogeneous 
stress 0, static electric as well as magnetic field or laser 
radiation (see [1] and references therein). The dynam- 
ics of dipolar excitons is complicated due to long-range, 
oc interaction in the system. In the present work, 

we demonstrate that under certain conditions the interac- 
tions in a dipolar exciton BEC can be considered as local 
interactions. As a result, the dynamics of the BEC can be 
effectively described by the generalized Gross-Pitaevskii 
equation (CPE). 

To characterize the non-stationary behavior of the 
dipolar exciton BEC we perform numerical simulations of 
the system where the resonant driving in the low- or high- 
frequency spectral domains is present. It was found that 
the dynamics of the system in these two cases are essen- 
tially different: for the low- frequency driving, the BEC 
spatial distribution tends to a stationary bell-like shape 
described well by the parabolic Thomas- Fermi profile [10| 
whereas for the high-frequency driving, the BEC demon- 
strates long-lasting, non-stationary oscillations. This os- 



cillatory state is somewhat similar to the wave turbulence 
state observed earlier in spatially restricted nonlinear su- 
perfluid systems H, 12| . 



Local approximation for a dipolar exciton BEC. At 
temperatures much lower than the BEC transition tem- 
perature, the dynamics of a dilute dipolar exciton con- 
densate is described by the generalized two-dimensional 
Gross-Pitaevskii equation 
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In Eq. ([T]) ^(r,t) is the exciton condensate wave func- 
tion, mex is the exciton mass, V{r) = is an exter- 
nal parabolic trapping potential, a is the trapping po- 
tential strength, U{p) = e^D^/ep^ is the pairwise ex- 
citon interaction potential, e is the electron charge, D 
is the inter-well distance, e is the dielectric permittiv- 
ity of the semiconductor, p = \r — is the distance 
between the excitons, and 7 = l/2rex, where Tgx is the 
exciton lifetime. Below an isotropic trap is considered. 
However, this is not a restriction of the model and we 
can also consider anisotropic traps with different trap- 
ping potential strengths in x and y directions. The last 
term in Eq. ([T]) describes creation of the excitons due 
to the interaction with the pumping laser radiation and 
the exciton decay similarly to that used in earlier works 
I3I . [1^ . However, to capture the resonant pumping in 
a given spectral domain, we introduce a linear operator 
R instead of the direct driving [13] or the homogenous, 
frequency-independent growth increment [14]. The R op- 
erator is defined via its matrix elements in a functional 
basis that enter into the system of equations (j4]) given 
below. 

The integral term in Eq. ([T]) diverges at small distances 
p, that corresponds to known divergence of the scattering 
amplitude of the dipolar excitons [4]. To regularize the 
integral we introduce a cutoff distance ro = {ep/ e^D'^Y^^ 
defined by the equation /7(ro) = /i, where p is the chem- 
ical potential in the exciton system. This is equiva- 
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FIG. 1: Steady- state spatial distribution for the dipolar ex- 
citon BEG density in a trap for the low-frequency driving at 
Ro — 0.1. Inset shows the evolution of the exciton density 
profile with rising the pumping rate Rq\ 0.1 (■), 0.3 (•), and 
0.5 (a). The curves represent the Thomas- Fermi distribution 
function. 



lent to the cutoff of the exciton energies at the chemi- 
cal potential ji. In this case, the function \^{r' ,t)\^ can 
be expanded in a Taylor series over p and the integral 
term reads j ^^^J^r'\^{r' ,t)\^U{p) = g\^{r,t)\^{l ^ 
O(ro/a)), where g = 27re^I^^/£ro, 0{ro/a) denotes the 
terms of (ro/a) order, and a ^ ^(r, t)/|V^(r, t)| is the 
characteristic length, at which the condensate wave func- 
tion changes significantly. For the low- frequency resonant 
driving a coincides by the order of magnitude with the 
size of the exciton cloud d 10 pm 0,[s|- Hence, the used 
approximation for the integral term is valid if ro <C 10 
jam. For the high-frequency driving, a can be smaller 
than d and is defined by the driving frequency. The es- 
timates for the experimental conditions 0, [sl show that 
ro/d < 0.04 and therefore, the corrections ^ (ro/a) are 
negligible. 

Thus, under the experimental conditions 0, Q the dy- 
namics of the dipolar exciton BEC is described by the 
equation 
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■A^(r,t) + l/(r)^(r,t) 



dt 2me 
+^^(r, |^(r, t)\^ ^ih(R- 7) ^(r, t). (2) 

Eq. (|2]) coincides with the "traditional" GPE, in which 
creation and decay of the particles are taken into account. 
However, in Eq. (|2]) the exciton interaction strength g 
depends on the chemical potential of the system and thus, 
on the density of the exciton BEC. In what follows we 
incorporate the equation for the chemical potential in the 
BEC determined for a dilute gas in a trap , 

p = +2i7(^))/iV, H^^^ = ^-j dr\^{r,t)\\ (3) 

where N = (1/2) / (ir|^(r,t)p is the total number of the 
excitons in the BEC, and and H'^^^ are the quadratic 
and fourth-order terms in the Hamiltonian of the system, 
respectively. 



Exciton BEC formation. To investigate the nonlinear 
dynamics of a dipolar exciton BEC, we solve the system 
(j2j), (j3j) with the help of the spectral representation that 
results in the following system of equations 

" -ig J2 Wr,,rr.,KlA*^AkAi + - 7)A„, (4) 



dt 



where Anit) are time-dependent spectral amplitudes, 
Wn^m,k,i is the matrix element of the interaction term 
H^^\ and Rn is the matrix element of the R operator. 
An(t) are the coefficients of the expansion of the con- 
densate wave function, ^(r,t) = ^71(^)^^71(^)7 where 
the basis functions ^nir) are the eigenfunctions of the 
Hamiltonian for a single particle in a parabolic poten- 



tial [l5|, and n = (nx^riy) is the two dimensional index. 



The Schrodinger equation for the basis functions coin- 
cides with the linear, Hermitian part (i.e. the first two 
terms in the r.h.s.) of Eq. (|2]). This choice allows us 
to capture the behavior of the condensate wave func- 
tion by taking into account a relatively small number 
of terms in the expansion for ^(r,t) and hence, obtain 
a good numerical convergence. In the simulations, the 
length and time are expressed in the oscillatory units 
£0 = {h^/arrie^Y^^ = 0.9 /im, to = (mex/o^)^/^ = 1.6 ns 
calculated for the external trapping potential at a = 50 
eV/cm^ and the exciton mass mex = 0.22mo, where 
mo is the free electron mass. The initial conditions 
were chosen in the form of quasi-equilibrium distribu- 
tion An{{)) = [T/ (/io + + ^2/ + 1)]"*^^^ exp(i^-„,) with the 
dimensionless temperature T = 0.1, the dimensionless 
chemical potential /io = 1? and random phases (p^. We 
found that the results only weakly depend on the choice 
of the T and /io constants. We numerically integrated 
the system of equations (|4|) with the 4*^ order Runge- 
Kutta scheme on a graphical processing unit NVIDIA 
Tesla S2050. 

We consider two cases: (a) the exciton system is 
pumped at low spectral modes, i.e., at frequencies com- 
parable with the fundamental frequency in the parabolic 
trap, ujQ = t^^] (b) the system is pumped in the high 
frequency spectral domain uo > uoq. It is worth noting 
that, in experiments the ratio of the driving frequency 
and the fundamental frequency of the trap may be tuned 
by changing both the laser pumping frequency and the 
trapping potential strength a 0,l8|. 

For the low- frequency driving, we set Rn = Rq n = 
{nl H- n^)^/^ < po where po is a given cut-off and Rn = 
otherwise. Below, we show the results obtained for 
Pq = 2. The driving was applied at t = 0. Initially, the 
system exhibits transient oscillations. At a later time, 
the oscillations were damped and the system relaxed to 
a stationary state, in which the spectral amplitudes An 
tend to constant values. The time required to approach 
the stationary state was equal to t ^ bOto ^ 80 ns for 
a low pumping rate Rq = 0.1 and it decreases to t ~ 
20^0 ~ 32 ns for a high pumping rate Rq = 0.6. We 
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FIG. 2: The exciton density profiles at t = 200to for (a) i?o = 0.1 and (b) Rq — 0.3 for tfie fiigfi- frequency driving at pi = 4, 
P2 = 6. Inset on frame (a) shows the exciton density plotted y — and averaged over the time period 50to < t < 200to and 
three independent runs (points). Curve shows the fitting by the Thomas- Fermi distribution. 



also observed that the initial exciton distribution relaxed 
to zero if the pumping rate was less that the threshold 
value Rq ~ 0.03. The presence of the finite excitation 
threshold is in agreement with the results of previous 
simulations of the polaritons dynamics [14]. Figure [1] 
shows the stationary non-uniform spatial distribution of 
the exciton BEC density at t = 200to for = 0.1. It 
is seen in Fig. [1] that an exciton "cloud" of size '^6 — 8 
£o that is, 5 — 7 fim is formed at the center of the 
trap. Formation of the cloud at the center of the trap is 
in agreement with the observations with indirect excitons 
in a trap [7] and in the BECs of trapped atomic gases [lo| . 
As is seen in the inset in Figdl the spatial distribution of 
the exciton in the BEC can by described by the Thomas- 
Fermi distribution = no[l — {r/d)'^] (no is the 
BEC density at the center of the trap and d is the effective 
radius of the cloud), which is used for characterization of 
inhomogeneous atomic BECs [10]. It is seen in Fig.[T]that 
the exciton density at the center of the trap gradually 
grows with the increase of the pumping rate. 
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FIG. 3: Dependence of the dipolar exciton interaction 
strength on time for the high-frequency driving at pi = 4, 
P2 = 6. Inset: Time dependence of the spectral amplitudes 
\An\ at n = (0,0) (the fundamental mode) and (8,8). 

When the system was driven in the high frequency 
spectral range {R = Rq dit pi < n < p2, and P2 > Pi > 4) 
the BEC exhibited persistent, non- damped oscillations. 
At relatively small pumping rates, Rq ~ 0.1, the spatial 
BEC profile at any given moment of time was far from 
the Thomas- Fermi distribution and it was represented as 
a series of irregular "spikes" , as seen in Fig. [2^. Never- 



theless, the mean exciton BEC density distribution aver- 
aged over a long enough period of time, (|^(r,t)p), was 
close to the Thomas-Fermi distribution, as shown in the 
inset in Fig. [2^. At larger pumping rates, Rq > 0.3, the 
density oscillated around a mean parabolic profile, which 
was close to the Thomas- Fermi distribution (Fig. ^jp). 
We emphasize that, unlike the case of the low frequency 
driving, the BEC oscillations were observed during the 
whole time period of the simulations 200to, which is much 
longer than the time duration of the initial transient pro- 
cesses. Figure [3] shows the respective oscillations of the 
spectral harmonics, A^,, at t > 50to and the variations 
of the dimensionless interaction strength, g' = gme^/h'^, 
with time. A decrease of both An and g' dX t < 50to 
corresponds to initial relaxation after the driving was ap- 
plied. 

To characterize the oscillatory regime, we plot a two- 
dimensional amplitude distribution (lA^P) averaged over 
time and three independent runs, as a function of 
{ux^riy), see inset in Fig. HI It is clearly seen that, de- 
spite the system is driven at the high-frequency range 
(marked by dashed curves), the maximum of the spec- 
tral amplitudes (deeper color) is positioned at low spec- 
tral numbers n < 2. Therefore, mutual interaction be- 
tween different spectral scales is present in the system. 
In other words, a flux of the excitations from the region, 
at which they are generated by an external pumping, to- 
ward the low-frequency region is formed. The resultant 
state is quite similar to a wave-turbulent state known for 
weakly nonlinear systems [13] • In those systems, non- 
linear interaction of running waves resulted in formation 
of turbulence, which was characterized by establishing 
of power-like, Kolmogorov spectra of energy and parti- 
cles distributions. In our case, the system is not spa- 
tially homogenous and the interacting, normal variables 
are present by the oscillatory modes An. 

To further characterize the turbulent regime, we calcu- 
late the angle- averaged distribution for the spectral am- 
plitudes, Nn, = Enltf^'d^^P)' ^^ig- SI Formation 
of power-like tails of the distribution Nn^ ex n'^ at fre- 
quencies below and above the characteristic driving fre- 
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FIG. 4: Angle- averaged occupation number of the excitonic 
modes in the BEG, Nn^, as a function of the radial number rir, 
plotted in log-log scale. The data are averaged over three in- 
dependent runs. The center of the pumping region is labeled 
by a vertical arrow, An = 3. The lines show a power- like dis- 
tribution for Nur- - const x at m = and m = —2. Inset 
shows the averaged spectral amplitude distribution |An,|^, the 
dashed curves show the boundaries of the pumping domain. 



quency is clearly seen on the main part of Fig. SI that 
supports our conclusion on the establishing of the turbu- 
lent regime in the dipolar exciton BEG in GQWs. The 
presence of two different exponents m = and m = — 2 
in the low- and high-frequency domains, respectively, in- 
dicate that the fluxes of the energy and of the number 
of particles through the scales are simultaneously formed 
in the system (cf. theory for traditional GPE jl6|). An 
essential difference between the previous considerations 
[lit results is that the effective damping in the 

system produced by the exciton relaxation is finite and 
cannot be disregarded for all modes. In effect, the fluxes 
of the energy and the number of particles are not com- 
pletely conserved in the spectral space. 

In conclusion, we demonstrate that the dynamics of 
the Bose-Einstein condensate of the excitons in cou- 
pled quantum wells can be described by the generalized 
Gross-Pitaevskii equation with the local (contact) inter- 
action, despite a long-range dipolar exciton interaction 
is present. The effective interaction strength ^ is a func- 
tion of the chemical potential of the system and therefore, 
should be determined self-consistently from the exciton 
distribution in the BEG. We show that, if the system is 
driven by an external pumping at low frequency modes, 
the spatial distribution of the excitons in a trap is de- 
scribed well by a parabolic density profile formerly known 
for dilute atomic BEGs. However, if the system is driven 
at high frequency modes, strong time-dependent density 
fluctuations are excited in the BEG, and the condensate 
density at any given moment of time can be far from 
the equilibrium parabolic profile. We infer that a tur- 
bulent state is formed in the exciton BEG in this case, 
which is characterized by a nonlocal particle balance in 
the system. The latter results in formation of the power- 
like spectra for the exciton distribution function. The 
turbulent state is somewhat similar to that recently ob 
served in superfluid ^He 



between all mentioned examples and the system under 
consideration is that in our case the interaction strength 
between the normal modes is a function of the occupa- 
tion number of the modes. Therefore, turbulence can 
not be considered as "weak" and it involves more com- 
plex interactions between the modes. It is worth not- 
ing that formation of long-living nonequilibrium states 
was recently described for a classical system with long- 
range interaction [l8[. The results of our consideration 
are useful for understanding of nonlinear phenomena in 
low-dimensional quantum systems including indirect ex- 
citons in GQWs 0,181, exciton polaritons propagation Q, 
and a photon BEG in dye microcavities [19|. 
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